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Physical  processes  in  blast,  perforation  and  penetration 

For  weapons  having  blast  characteristics  alone,  the  loading  usually 
engulfs  the  strut'.'. are  and  can  produce  lethal  damage  through  overall 
deformation  of  the  vehicle  or  through  damaging  of  a  critical  section- 
For  fragment  typ«  weapons  the  high  velocity  fragments  perforate  parts 
of  the  structure  and  can  produce  lethal  damage  if  critical  portions 
of  the  structure  are  blown  out  by  the  fragments.  There  are  other 
weap. is  which  have  both  fragment  and  blast  characteristics  and  these 
can  produce  structural  damage  by  the  fragments  first  taking  out  part 
of  the  structure  or  plastically  deforming  the  structure  (without  per¬ 
foration)  and  then  the  blast  finishing  the  damage. 

A  description  of  the  processes  involved  in  penetration  and  perforation 
of  very  high  velocity  fragments  against  plate  type  targets  is  contained 
in  an  earlier  reference.1  When  a  high  speed  fragment  or  projectile 
reaches  the  surface  of  a  target  several  things  can  occur  depending  upon 
the  characteristics  of  the  fragment  and  target.  The  fragment  can  per¬ 
forate  the  target  completely  and  continue  traveling  with  a  small  change 
in  velocity  due  to  the  resistance  of  the  target.  If  the  target  is  an 
aircraft  or  missile  structure  the  fragment  can  enter  on  one  side  and, 
if  the  speed  is  high  enough,  and  the  target  resistance  low  enough,  it 
will  leave  on  the  other  side.  This  will  happen  in  very  thin  structures 
which  offer  little  resistance  to  perforation  such  as  shown  in  Figure  1. 
For  thicker  targets  offering  more  resistance,  the  fragment  could  per¬ 
forate  the  entrance  side  and  experience  a  substantial  decrease  in  velo¬ 
city  —  thus  having  its  kinetic  energy  considerably  lowered.  Part  of 
this  energy  change  will  go_intib  perforating  and  deforming  the  target 
and  projectile  anti  part  will  go  into  heating  the  fragment  and  target 
possibly  to  the.  point  of  vaporization.  In  the  meantime,  as  perforation 
occurs,  liquid  spall  fragments  can  he  formed  on  the  inside  surface. 

This  physical  process  is  illustrated,  in  Figure  2.  In  some  cases  the 
entering  fragment  will  not  be  vaporized  because  the  initial  velocity 
is  not  large  enough! . In  this  case  the  spall  fragments  together  with 
the  original  broken  up  fragment  will  load  the  exit  side  as  shown  in 
Figure  3.  If  the  velocity  of  the  original  fragment  is  not  high  enough 
for  perforation  then  penetration  could  occur  as  shown  in  Figure  4.  If 
the  velocity  is  not  sufficient  for  penetration,  then  the  fragment  might 
cause  pla-tic  deformation  of  the  structure  and  bounce  off  such  as  shown 
in  Figure  5. 

In  the  process  illustrated  in  Figure  2  the  shock  will  undoubtedly  blow 
out  the  entrance  side  of  the  structure  since  the  energy  of  the  enter¬ 
ing  figment  has  been  converted  to  heat  during  the  perforation  process 
and  the  fragment  will  he  converted  mostly  to  liquid  and  gas  at  the  rear 
of  the  entrance  side.  In  the  process  shown  in  Figure  3,  the  spall 
fragments  will  probably  blow  out  the  exit  side  unless  their  energy  is 
absorbed  in  internal  components.  In  this  case  fragments  v/ill  be  formed 
at  the  inside  of  the  entrance  panel  and  be  projected  toward  the  e’-tit 

♦Superscripts  refer  to  references  listed  at  the  end  of  the  report. 
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Fig.  1  Perforation 


Fig.  2  Melting  and  Vaporization 


Fig.  3  Breakup 


Fig.  5  Plastic  Deformation  and 
Bounce  Off 
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side.  Liquification  and  vaporization  will  ordinarily  occur  at  velo¬ 
cities  of  impact  greater  than  25,000  ft. /sec.  for  steel  fragments  per¬ 
forating  an  aluminum  target.  For  any  velocities  lower  than  this,  the 
phenomena  described  in  Figures  2.3,4  or  5  will  probably  occur. 


II.  Energy  concepts  as  a  general  description  of  damage 

In  blast  work  one  of  the  most  com  on  ways  to  show  damage  results  is 
on  an  isodamage  curve.  This  iso  carnage  curve  consists  of  a  pressure- 
impulse  plot  such  as  shown  in  Figure  6 


,  ? 

It  can  be  shown  theoretically  that  the  curves  of  constant  damage 
plotted  in  this  P-I  plane  are  approximately  constant  blast  energy 
curves  since  the  blast  energy  flux  (energy  per  unit  area)  can  be 
written 


where  E  -  energy  flux  in  the  explosion  at  the  target 
B 


I  =  Impulse  per  unit  area  in  the  blast 


[1] 


P  =  blast  overpressure  at  the  target 
A-  **  density 
C0=  sound  velocity  in  the  air 


If  the  energy  available  to  do  damage  (i.e.  energy  f ±,a)  is  equated 
to  the  energy  absorbed  in  the  structure  then  the  curve  plotted  in 
the  P-I  plane  is  both  a  constant  damage  curve  as  well  as  a  curve  of 
constant  available  blast  energy. 

If  we  now  consider  a  weapon  in  which  fragments  are  a  mechanism  of  dam¬ 
age,  the  available  energy  from  each  fragment  is 

E  -  £  [2] 

where  M  =  mass  of  fragment 

rT  -  velocity  of  fragment  at  target 


A  curve  plotted  in  the  v  -  K  plane  will  look  like  the  curve  shown  in 
Figure  7 
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Fig.  7  Isodamage  Fragment  Curve 

By  the  same  reasoning  used  for  blast,  the  curves  plotted  in  the 
v  -  M  plane  for  fragments  are  also  constant  energy  curves  and  the 
isodamage  curve  for  fragments  is  the  v  -  M  curve. 


III.  Basic  Equations  for  Fragment  Perforation  and  Collapse  of  Cylindrical 
Structures 


A.  Fragmentation  -Perforation 

If  the  mass  and  velocity  of  the  incoming  fragments  are  of  suffic¬ 
ient  magnitude  then  these  fragments  will  perforate  the  structure. 
The  Thor  equation  for  determining  the  change  in  velocity  for  frag¬ 
ments  which  perforate  a  target  is  as  follows: 

Vh  -  VS  -  /oe  [3] 

where 

Vr  is  the  fragment  residual  velocity  in  fps  after  perforation 
has  occurred 

Vs  is  the  fragment  striking  velocity  in  fps 
e  is  the  target  thickness  in  inches 

A  is  the  average  impact  area  of  the  fragment  in  square  inches 

is  the  weight  of  the  original  fragment  in  grains 

6  is  the  angle  between  the  trajectory  of  the  fragment  and  the 
normal  to  the  target  material 

C.otfi  <r  A  are  constants  determined  for  each  material  and  are 
shown  in  Table  Is 
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Table  1  Perforation 


Residual  Velocity  Constants 

(Steel  Fragments) 

Material 

C 

•<  -/ 

> 

Magnesium 

6.3 

1.1  1.2 

1.1 

.09 

Aluminum  Alloy 

7.0 

1.0  1.1 

1.2 

.14 

Titanium  Alloy 

6.3 

1.1  1.1 

1.4 

.17 

Cast  Iron 

4.8 

1.0  1 ..  1 

1.0 

.52 

Face  Hardened  Steel 

4.4 

.7  .8 

1.0 

.43 

Mild  Steel 

6.4 

.9  .9 

1.3 

.02 

Hard  Steel 

6.5 

.9  .9 

1.3 

.02 

Copoer 

2.8 

.7  .7 

.8 

.80 

Lead 

2.0 

.5  .5 

.7 

.82 

Tuballoy 

2.5 

.6  .6 

.9 

.83 

For  fragments  which 

are  on 

the  threshold  of  perforation  (i. 

. e.  if  the 

frcgment  had  a  lower 

velocity,  it  would  not  perforate) 

the 

velocity 

given  by 

Va  =  1C 

.  ✓  .  \*^l  fil  s-  „  \ 

1  (e/U  ryts  ^^=9) 

where  the  constants 

are  given  in  Table  2. 

Table  2  Perforation 

Threshhold  Velocity  Constants 

(Steel  Fragments) 

Material 

*/ 

-A 

K 

Magnesium 

6.4 

1.0  1.1 

1.0 

Aluminum  Alloy 

6.2 

.9  .9 

1.1 

Titanium  Alloy 

7.6 

1.3  1.3 

1.6 

Cast  iron 

10.2 

2.2  2.2 

2.2 

Face  Hardened  Steel 

7.7 

1.2  1.4 

1.7 

Mild  Steel 

6.5 

.5  1.0 

1.3 

Hard  Steel 

6.6 

.9  1.0 

1.3 

Copper 

14.1 

3.5  3.7 

4.3 

Lead 

11.0 

2.7  2.7 

3.6 

Tuballoy 

14.8 

3.4  3.5 

5.0 

A  fragment  with  a  velocity  slightly  greater  than  this  velocity  (eq.  [4]) 
will  knock  out  a  section  of  the  structure  of  thickness  e  and  area  A. 

For  a  given  fragment  distribution  in  which  we  have  effective  area  A^ 
over  which  fragments  perforate  the  weapon  will  knock  out  this  area  of 
the  structure. 

Equations  [3]  or  [4]  are  used  only  in  one  way.  Given  the  mass  of  the 
fragment  ms,  its  projected  area  a,  the  target  thickness  e,  the  angle  e, 
and  the  appropriate  constants,  this  equation  gives  the  minimum  velocity 
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of  fragment  which  will  cause  perforation.  Let  the  velocity  of  the 
fragment  be  V.  If  V<V0  no  perforation  will  occur,  if  V  >  VQ  perfora- 
tion  will  occur  and  the  fragmentation  damage  will  be  done.  If  the  mass 
and  velocity  parameters  of  the  warhead  and  the  thickness  of  the  target 
are  of  appropriate  value  so  as  to  cause  perforation  over  a  large  enough 
area,  then  perforation  could  be  sufficient  to  kill  the  structure  alone 
without  any  further  considerations. 


B.  Fragments  -  No  perforation 


C. 


If  fragments  do  not  perforate  they  still  could  inflict  some  damage  on 
the  structure  by  imparting  their  energy  to  the  structure  and  causing 
some  plastic  deformation  by  penetrating  the  skin  and/or  hitting  the 
structure  and  bouncing  off.  Let  the  kinetic  energy  of  a  fragment  (which 
either  penetrates  or  bounces  off)  of  mass  ms  and  velocity  /vf be  denoted 
by  E3,  then 


■  i  [5] 

Assuming  no  energy  Josses,  the  total  energy  imparted  by  these  fragments 
will  be  the  sum  over  s,  i.e.* 


.S' 


s  < 


[6] 


Collapse  in  a  cylindrical  structure  -  Missile  bodies  and  fuselages 


1.  Collapse  due  to  non  perforating  fragments 


It  is  conceivable  that  the  totai  energy  from  non  perforating  frag¬ 
ments  could  be  sufficient  to  cause  collapse  in  a  cylindrical  shell 
structure  such  as  shown  in  Fig.  8  and  9  **  Assuming  a  linear  har¬ 
dening  material  such  as  shown  in  Figure 10,  the  energy  absorbed  in 
the  collapse  mode  can  be  written:^ 


\7=  j/T  A/ 


A  j  r,  +X  v,  -  A/T  7rCs  J 


\  3 


i 


C  t-J' 


where  =  yield  stress  of  material  in  pure  tension 
h  =  thickness  of  material 
1  =  length  of  shell 
A  =  hardening  parameter 
es  =  yield  strain 
■J  -  Poisson's  ratio 


[7] 


Temperature_  effects  are  considered  by  adjusting  OJ  and  X  .  The 
values  of  X,  }  v,  are  given  in  Fig. 11  The  above  relation,  [7]  is 
the  energy  necessary  to  cause  collapse  with  ?.  given  deflection 
as  shown  in  Figures  9  and  11 

*In  order  to  assume  that  the  energy  from  the  fragments  is  imparted  to 
the  structure  when  there  is  no  perforation, ' the  value  of ^  must  ap¬ 
proach  Vo  (the  perforation  velocity) ;  otherwise  the  fragments  will 
just  bounce  off  and  only  impart  a  small  portion  of  their  energy 
thereby  causing  only  elastic  deformation  in  the  structure. 

**Since  we  are  primarily  concerned  with  focused  fragment  -  blast  war¬ 
heads  in  this  report,  the  nature  of  the  loading  will  eliminate  the 
lobe  buckling  type  of  failure  which  is  characteristic  of  blast  load 
engulfment. 


11 

}j 


2.  Collapse  with  non  perforating  fragments  together  with  blast 

It  is  also  possible  that  a  combination  of  non  perforating  fragments 
and  blast  could  cause  considerable  damage  even  if  no  perforation 
occurs.  This  is  accomplished  by  energy  transfer.  First  the  frag¬ 
ments  impart  their  energy  E^  and  deform  the  shell  plastically  as  de¬ 
scribed  above.  The  blast  then  comes  along  and  imparts  its  energy  Efi. 

For  a  given  overpressure  P,.  the  value  of  bjf4C0  is  obtained  from 

T  Ps*- 

Figure  12.  E^  is  the  energy  per  unit  area  in  the  blast,  is  the 
air  density,  CD  is  the  sound  velocity  in  the  ambient  air,  T  is  the 
positive  duration  of  the  overpressure,  thus 

where  A  is  the  area  of  the  target  facing  the  blast.  The  total  ener¬ 
gy  imparted  to  the  target  by  non  perforating  fragments  plus  the  blast 

Et.Ui  ^  Ee  _  t91 

It  is  this  energy  (  )  which  is  compared  with  V  to  determine 

what  degree  of  collapse  (i.e.  what  )  the  combined  weapon  will 

inflict  on  the  target. 

3.  Collapse  with  perforating  fragments  and  blast 

If  the  fragments  have  already  perforated  then  they  have  knocked  out 
a  given  area  of  the  shell.  This  means  that  the  rest  of  the  shell 
is  only  capable  of  absorbing  energy  over  its  remaining  area.  The 
energy  necesasary  to  inflict  collapse  damage  with  deflection 
on  the  perforated  structure  is 

-  _  -  Eemaining  area  of  the  body  after  Penetration  takes  place  .  _ 
VP  ~  Original  area  of  the  body 


We  compare  Vp  with  Efi  (the  energy  available  frcm  the  blast)  to  see 

how  much  more  damage  will  occur  in  the  perforated  shell. 

4.  The  combined  case  -  see  Figure  13 

For  most  cases  some  of  the  fragments  will  perforate  thus  weakening 
the  structure,  the  blast  will  then  come  along  and  further  deform 
the  structure  and  then  the  slower  fragments  will  either  penetrate 
or  bounce  off  and  plastically  deform  the  structure.  The  energy 
available  from  the  weapon  to  damage  the  perforated  structure  is 


E /*,  •f-  s  J  ««  r  fr  * 


This  is  the  energy  to  be  compared  with  Vp  in  order  to  determine  how 
much  damage  this  available  energy  could  inflict  on  the  perforated 
structure. 


-7 


Schematic  of  Various  Types  of  Failure 
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Fig.  9  •  Collapse  Pattern 
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Fig.  9a  Buckling  Pat. tern 
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Figure.  10  Elastic-Linear  Hardening  Law 
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After  Perforation 


Residual  Energy  Available  from  Weapon 
(non-perforating  frags  +blast) 

Energy  Absorbed  in  Weakened  (Perforated) 
Section 


P*  netrated 


Perforated 

/  / 
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Pig.  13  Combined  Damage 


IV.  Buckling  Failure  of  Wings  and  Fuselages 


It  is  assumed  that  failure  both  in  a  wing  (or  other  control  surface) 
and  in  a  fuselage  takes  place  by  local  buckling  of  the  skin,  stringer, 
and  spars  such  as  shown  in  Figures  14  -  16  The  section  at  which  failure 
will  o^cur  will  be  the  section  which  is  weakened  by  fragment  perfora¬ 
tion.  In  order  to  determine  whether  steel  fragments  will  perforate 
aluminum  structures  we  merely  have  to  apply  the  THOR  equations.  After 
perforation  takes  place  the  local  section  will  be  weakened  because 
it  has  less  load  carrying  area.  This  remaining  -weakened  structure 
will  then  be  exposed  to  the  oncoming  blast. 


The  structural  theory  of  Rost  fui lure  behavior  of  buckled  structures 
was  worked  out  by  D'Amato'>  some  years  ago.  Neglecting  any  elastic 
vibrational  energy,  the  energy  absorbed  by  the  structure  is  the  elas¬ 
tic  energy  in  taking  the  structure  up  to  failure  at  the  weakened  cri¬ 
tical  section  plus  the  plastic  energy  in  rotating  the  structure 
through  angle  6  (see  Figurel7)  around  the  critical  section. 


Fig.  17  Configurations  of  the  Structure 

At  the  critical  section  there  is  a  resisting  moment  which  is  schematic¬ 
ally  shown  in  Figure  18 
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Figure  18  Hinge  Resisting  Moment 

The  area  under  the  M  -  &  curve  is  the  strain  energy  absorbed  as  the 
structure  undergoes  failure. 

Considering  the  structure  as  a  single  degree  of  freedom  deforming  in 
the  hinge  mode  (see  Fig.  4),  the  equation  of  motion  is  as  follows: 

TO  -  Mt'&)  -  M(&)  [U] 

where  0  =  rotation  angle 

I  -  mass  moment  of  inertia  of  the  hinged  section  around  an 
axis  through  the  critical  section 

,  ('■/■)  =.-  external  moment  applied  to  the  weakened  (perforated)  struc¬ 
ture  by  the  blast 

M  f<9)  =  resisting  moment  at  the  critical  section 

Following  the  logic  of  Ref.  2,  the  external  work  done  by  loading  up 
to  time  t  is 

.  eft)  *  t 

w(t)  -  J  - jtl  (i )  yn  <Jt  [12] 

from  [ID  ,.h  ° 

J-  f  r  .. 

iit  x  J  [n(t)-nre)ldt  U3] 


Thus 


Wd)  -  Ja  j  [M  ft)  -  Mrs)]  aft}  cH- 


A 


ill 


m 


In  order  to  calculate  the  total  work  done  this  .last  expression  has  to 
be  integrated  up  to  che  time  at  which  the  maximum  deflection  takes  place. 
If  the  time  at  which  this  maximum  deflection  taV.es  place  is  much  grea¬ 
ter  than  the  time  of  duration  of  the  load,  T,  then  we  o^y  have  to  in¬ 
tegrate  up  to  T.  Under  such  circumstances,  rtf&)  will  be  small  com¬ 
pared  to  AW  i  during  the  time  0  to  T.  Practical  structures  such  as 


[15] 


fuselages  and  wings  under  blast  at  distances  of  10* -  50'*  will  fall 
into  this  catagory.  Thu3 

W  =  f  V?  Ct  )j  j:  J  dt 

or 


where 


H  x 

W  -  iX 

ZT 


[16] 


H  -  /  MC+)Jt  fl7J 

Jo  L  ^ 

The  energy  available  from  the  blast  field  is  ™  /z  I  (i.e.  the  external 
work)  and  this  is  equated  to  the  energy  absorbed  by  the  structure, 

which  is  the  internal  work. 

For  a  flat  section  such  as  a  wing 


CVi'/zf*/  Scc’f'c** 


Prc  t  Cc^itr 
T 

H  - f t pet)  A  Celt  -  Ac  T 


[18] 


where 


and 


pCt  ) 
A 
C 
2 


=  blast  pressure  (assumed  uniform) 

=  surface  area  of  wing  from  critical  section  to  edge 

*  distance  from  critical  section  to  pressure  center 

*  f  -pCt)  &i  t  -  applied  impulse 

n  A 


3 


Thus  for  the  wing  or  other  flat  control  surface 

*($)'(&+ ) 


[19] 


[20] 


For  the  fuselage,  the  loading  tending  to  produce  moment  is  acting  only 
over  half  the  surface  area  (i.e.  over  the  half  facing  the  blast).  Thus 

H  ~  p(t)  ^  --^C  -T  [21] 


±C- 

zi  " 


3 

1 


<$)'(-£*) 


[22] 


In  cases  where  the  wing  or  fuselage  is  loaded  with  a  static  loading 
such  as  fuel,  end  tanks,  etc.  there  is  an  additional  turning  moment  pre 
sent.  Once  the  structure  is  weakened  by  the  fragments  taking  out  mater 
ial  at  a  given  section,  the  blast  arrives  and  can  do  additional  damage 
(i.e.  produce  a  larger  hinge  angle  around  the  critical  section)  if  it 
contains  sufficient  energy.  After  this  deformat  ton  is  completed,  the 
static  moment  can  then  produce  further  hinging  if  it  is  of  sufficient 
magnitude.  This  type  of  loading  will  be  discussed  in  the  sample  pro¬ 
blem  given  later  ^n  this  report. 


*i.e.  for  f rag-blast  type  weapons 
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V  .  Calculation  of  the  internal  work 
A.  Flat  surfaces 

The  internal  work  consists  of  the  elastic  work  done  on  the  struc¬ 
ture  bringing  it  up  to  failure  plus  the  plastic  work  done  by  the 
hinqe  restraining  moment  at  the  critical  section.  In  wings,  ele¬ 
vators  and  other  control  surfaces  which  are  cantilevered  from  the 
fuselage,  the  failure  takes  place  by  way  of  a  hinging  at  a  weaken¬ 
ed  section.5  This  weakened  auction  is  formed  when  a  portion  of 
the  structure  such  as  skin  or  spar  caps,  etc.  is  knocked  out  by 
fragments.  If  the  fragments  perforate  across  the  entire  surface, 
thus  cutting  the  surface  into  two  parts,  then  this  will  consti¬ 
tute  kill  by  itself  without  any  further  calculations.  If  only 
a  portion  of  the  section  is  knocked  out  we  can  then  calculate 
how  much  additional  damage  will  be  done  by  blast. 

Given  a  typical  wing  section  such  as  shown  in  Figure  19 


Fig.  19  Typical  Wing  Section 


First  the  failure  moment  is  computed  in  accordance  with  the  re¬ 
lation 


n 


[23] 


where  I  is  the  area  moment  of  inertia  of  the  weakened  section 
(a ‘ter  perforation  takes  place)  and  h/2  is  the  half  depth  of  the 
control  surface.  /  7_  , 

£  ft* 


dV-r  A. 


icLB.  {&  V 

/  o  * )  (  hs  J 


[241 
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rfhere  is  the  plasticity  reduction  factor  obtained  from  a  previous  re¬ 
ference,^  E  is  the  modulus  of  elasticity,  -*)  is  Poisson's  ratio,  tg 
is  the  thickness  of  the  cover  skin  (see  Fig. 19),  bs  is  the  width  of 
the  cover  skin  between  rivet  lines  attaching  skin  to  flanges  and 


where  bw  =  h 

-  web  thickness  (see  Fig. 19) 

f  =  distance  from  centerline  of  web  to  ,he  point  where 
tb'  'et  clamps  the  flange  to  the  cover  skin 


After  failure  has  occurred,  the  moment  -  hinge  angle  relation  is  given 
by  the  following  equation: 

+fi)  +Crviot  ~  71  fa  [2b] 

where  ?Y)ew  -  plastic  moment  per  unit  length  of  the  web  ~ 


in  which  yield  stress  of  web 

=  thickness  of  web 
h  s®  depth  of  beam 

n  =  number  of  webs  in  cross  section 
c  «  width  of  cover  skin 


7*1*5  =  plastic  moment  per  unit  width  of  cover  skin  = 
where  0^s  -  yield  stress  of  cover  skin 

ts  **  thickness  of  cover  skin 

5 

r  =  web  folding  radius  (see  D' Amato  ) 

The  functions  \fr^  y“'f '  are  dependent  upon  the  hinge  angle  9 

and  are  as  follows: 


(<f>  \  /  )"* 

<P C /  t-c* «&  +  f<u~. 

&  f  c<p  *  -  ^  a  )(t  -  & ) 


/  ^  Ccj  & 


i\/e) 


t  )  w  &  /  2  cf  c^oX  {9 

•  { <f  z-  0  &  ~  2  cf}  J  t  *  V.  c9 


<?,  ( I4 

/  -  f  **■»  - 


■dS.  <%  j 
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£&*■  zdr 

St 


-~~±^y~SL.  ( ep  c^  &/. 

S  4  4-*+\\.  Si 


4  *^  4i } 


cf  -  S/*A 


Si  =■ 
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_ /_ _ 

^  /•  f«a  <5^ 
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The  elastic  energy  absorbed  by  the 


structure  can  be  written” 


2.  El 


[28] 


where  M  is  the  bending  moment  in  the  stsucture, 
elasticity  and  I  is  the  area  moment  of  inertia. 


% 

A 


UuJl 


ft 

n  k 

zi 

H 


E  is  the  modulus  of 


For  a  uniformly  loaded  cantilever  beam  such  as  shown  above,  the  bending 
moment  is  as  follows: 

n=h 

The  maximum  moment  will  occur  at  the  root  and  will  be 


-  4  ft  f*" 

Thus  .  /,-  Nx. 

n  -  m.  (x/j) 

[29] 

and  the  energy  absorbed  elastically  will  be 

[30] 

The  plastic  energy  absorbed  will  be  the  area  under  the  M  -  9  curve  in 
the  plastic  region  shown  in  Figure  18 


B.  Cylindrical  shells 

5 

We  can  use  D1 Amato's  theory  for  the  fuselage  by  making  an  equivalent 
flat  structure  from  the  fuselage.  To  do  this  we  equate  the  load  carry¬ 
ing  area  in  the  top  and  bottom  plates  of  a  flat  structure  to  the  load 
carrying  area  of  the  cylindrical  shell  of  the  same  thickness. 
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Equivalent  Flat  Structure 

Thu.  Zt*r**Krt  (31) 

-*cr  ~  7 rr 

The  width  of  the  flat  structure  will  be  h  of  the  circumference  of  the 
cylinder. 


The  height, h, of  the  equivalent  flat  structure  is  found  by  equating 
the  area  moments  of  inertia  of  the  f1at  structure  and  cylinder 

thus  2  '  ^2.  J  '  r^r  2  ^ 

and  .  ,  [321 

K  -  {J  2.  r 

The  equivalent  thickness  of  the  fuselage  which  contains  plating  and 
attached  stringers  will  be  taken  as 

i  -  AnL  J-  t 

* c  ~  ~rro  °  {33] 

where  A  *  cross  sectional  area  of  stiffeners 

n  “  number  of  longitudinals  in  cross  section 
te  -  euqivalent  uniform  thickness  of  cylinder 
t0  *  thickness  of  skin  itself 

VI .  Sample  problem 

A  typical  wing  and  fuselage  section  were  considered.  The  blast  energy 
available  for  6.5#  of  Explosive  as  a  function  of  distance  is  aiven  in 
Fig. 20  for  the  wing  section.  Impulse  values  given  by  Goodman4  were 
used.  The  energy  absorption  (or  internal  work)  curves  for  the  struc¬ 
ture  are  shown  in  Figure  21  .  In  this  particular  case  there  was  a  static 
loading  on  the  wing  having  a  moment  of  1,000,000  #  in.  around  the  cri¬ 
tical  section  (Sta.  143).  Figs.  20  and  21  show  that  the  energy  content 

of  the  6.5#  blast  is  so  small  that  it  would  do  practically  no  damage 
to  the  wing  even  if  8096  of  the  wing  was  knocked  out  by  fragments.  Fig¬ 
ures  22  and  23  show  that  the  static  moment  is  so  large  that  it  will  fail 
the  wing  after  the  fragments  have  done  their  damage..  The  failure  moment 
is  shown  in  Fig. 22.  This  figure  shows  that  the  static  moment  will  not 

inflict  damage  until  about  70"  of  structure  is  knocked  out  at  the  cri¬ 

tical  section.  Fig.  23  illustrates  that  the  failure  will  be  complete 
since  the  resisting  moment  for  a  structure  with  70"  removed  is  less 
than  half  of  the  static  moment. 

Fig.  24  illustrates  the  energy  available  in  the  6.5#  blast  on  the  fuse¬ 
lage  and  Fig.  25  shows  the  internal  work  curves  for  t.he  fuselage  section 
near  the  tail.  These  curves  show  that  if  most  of  the  circumference  is 
knocked  out  by  the  fragments  then  the  6,5#  blast  very  close  to  the  fuse¬ 
lage  (less  than  5  feet)  might  produce  a  small  post  failure  hinge  angle. 
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trii.  Determination  of  modes  of  failure  for  monocoque  structures  under 

pure  blast 

Experiments  have  shown ^  ^  that  monocoque  cylindrical  shells  sub¬ 
ject  to  side  on  blast  will  fail  in  one  of  two  types  of  modes.  The 
first  is  the  collapse  mode  as  seen'1  in  Figure  9  and  the  second  is 
lobar  buckling  as  seen  in  Figure  9a,  There  is  a  controversy  as  to 
whether  the  collapse  is  an  instability  known  as  collapse  hinge  buck- 
3ing-*"}r'  ,  or  whether  it  is  a  straight  collapse  determined  by  the 

criterion  that  the  Von  Mises  yield  condition  is  satisfied  at  the  max¬ 
imum  stress  pints.  Figure  26  shows  curves  computed  for  the  three 
types  of  failure.  The  curves  for  collapse  hinge  buckling  and  lobar 
buckling  were  taken  from  Reference  9.  In  these  curves  L  is  the  shell 
length,  a  is  the  radius  and  h  is  the  thickness.  The  collapse  curves 
were  computed  using  membrane  theory2  as  follows: 

If  f\f<p  is  tht  peripheral  membrane  stress  resultant.  A/*  the  longi¬ 
tudinal  membrane  stress  resultant  and  A/*#  the  membrane  shear 
stress  resultant  then 


-  a. 


■p  (<p  ) 


where  is  the  peripheral  pressure  distribution 
N*  »  -i£  0?  u-  4-X.  V  ^ 
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Let  the  length  of  the  shell  be  L  and  assume  that  it  is  supported 
by  a  diaphragm  at  each  end  X.  -  +  L /z  .  if  the  origin  is 

at  the  center  of.  the  shell  the  boundary  conditions  are 

Nx.  -  O  <*>  X.  -  ±  >-/ 2. 


[35] 


The  pressure  used  in  computing  the  curves  of  Ref.  9  (Fig.  26) 
was 


sp  C<?)  -  ( >°K  ~Px)c<^cf  *  PZ  -J 


=  Px 


E<<$*  Irr 


[36] 
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Thus 


A/<^>  =  G~CPft-~?r  V  ^  <2-  Pr 

Ffcp)  -  -  (Pit  -  Px  Zcp  f  3  7  J 

A/x^>  “  X  (Pit-  Px  )++**■  ?-'P 
Nx.  -  ^  (f  ^  *)(?*  -  Pt)Z  Ccr^  2  cf 

Plasticity  will  first  begin  at  X  -<P-o  if  the  Von  Mises  yield  condition 
is  satisfied  at  this  point,  i.e.  ' 

A /j  1  —  A/*  Ncfi  +  fV<p  ~h  3  /V xtf  ~  -A-  t  J 

where  (7J  =  yield  stress  in  pure  tension 
Thus  the  equation  for  the  determination  of  the  critical  value  of  P/tCc'Pz) 


F^P'.-pxf-Fp^-p*)  +  p^-- 

Sc4  p  4*  a 

Given  the  ratio  '^/Px  ,  the  critical  value  of  Px  can  be  determined.  The 
value  of  q  plotted  in  Fig.  26  is 

t "It"1'  [4C 

The  values  computed  from  the  above  equation  for  Pt./pT  :  2^  T  are  shown 
in  Figure  26  along  with  the  values  for  collapse  hinge  buckling  and  lobar 
buckling  given  in  Reference  9.  The  co? lapse  curves  were  computed  using 
,)  =  .3  and  ^/g~s  -  1000,  which  seemed  realistic  for  steel  cylinders. 

Among  Schuman's  tests^-®?^  there  were  seven  steel  cylinders,  each  with 


=  79  and  ^Pol  varying  from  4  to  16.  All  seven  cylinders  went  into  the 
collapse  mode  of  failure.  Examination  of  Figure  26  reveals  that  for  L3x. 
between  4  and  16  collapse  would  occur  before  lobar  buckling  and  lobar 
buckling  would  occur  before  collapse  hinge  buckling.  Since  the  tests 
showed  collapse  patterns,  it  can  be  concluded  that  all  the  cylinders 
failed  in  ordinary  collapse  as  described  by  the  theory  in  this  section 
of  the  report.  Extrapolation  of  the  curves  for  aPP  =79  from  the  value 
of  ^/PL  =  64.6  shown  in  Figure  26  will  change  the  curves  only  slightly. 
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VIII.  Critical  non-dimensional  parameters  in  fragment-blast  response 
A'.  Perforation 


In  section  III  A  the  basic  empirical  equations  for  perforation 
were  given.  A  theoretical  relation  for  residual  velocity  is 
given  in  Goldsmith's  book^"*  and  is  as  follows; 

Vr  -  vs  e 


[41] 


where 


residual  velocity  after  penetration 
initial  velocity 

mass  density  of  plate  being  perforated 
frontal  area  being  perforated 
thickness  of  plate  being  perforated 


'yn  =  mass  of  the  projectile  doing  the  perforating 

The  quantity  /  A  A0  is  the  mass  of  the  plug  which  is  taken 
cut  of  the  plate.  Therefore  the  above  equation  becomes 


Vr  &  ~ 
V,  "  L 


[42] 


Thus  the  ratio  of  residual  to  initial  velocity  is  dependent  upon 
the  ratio  of  the  mass  of  the  plug  taken  out  of  the  plate  to  the 
mass  of  the  projectile  perforating  the  plate. 

Figure  27  shows  a  comparison  between  the  results  obtained  from 
equations  [3]  and  [41 j  .  it  is  seen  that  the  results  compare 
reasonably  well.  Thus  the  important  dimensionless  ratio  in  de¬ 
termining  residual  velocity  is  . 


Fragments  Penetrating  Aluminum  Plate 
Initial  Velocity  =  5800  ft. /sec. 

Mass  of  Projectile  =  31  grains  (V  cube) 


Plate- inches 


B.  Post  failure  buckling 
1.  Fuselages 

In  section  VB  of  this  report  we  assumed  that  once  the  weakened 
section  was  formed  by  perforation,  the  fuselage  would  then  buckle 
at  the  weakened  section.  The  moment-hinge  angle  (M-e) relation  fcr 
the  fuselage  was  determined  by  the  equation 


(from  eq.  [ 26] ) 


[43] 


w  v  C' fa 

where  C  7Tr  (from  eq.  [31]  ) 

^  ^sl/4 

^  -  A  r  .  4 

ts  * 

From  equation  [27]  it  is  seen  that  the  value  of  )^j  depends  upon 
the  hinge  angle  6  and  cfi  where  A*  is  the  wave  length 

of  the  buckle  and  h  is  the  depth  of  the  beam.  The  value  of  ^ 
for  various  values  of  F  (=  (p  )  is  given  in^Table  2  along  with  the 
values  of  ^3  integrated  over  &  (i.e.  f  )J &  ). 

The  energy  absorbed  at  the  hinge  is  * 

E  =  £Vwie j  J  yr3  CFJd)d&  [44] 

*  Cm 


Using  Table  2  and  the  relations  given  previously,  the  M-e  curve 
and  the  energy  absorption  can  be  readily  computed. 


2.  Wings 

For  the  wing  section  the  post  failure  moment  depends  upon  \jf 
)  and  the  energy  absorption  at  the  hinge  will  be 

£  -  +  +  ^  t45J 

Table  2  contains  1 \Jr  J jr  yr  as  functions  of  F(=^  and  0  as 
well  as  /yritlQ  j8  ,L.  ^  Cindd  •  The  nomenclature  in  Table  2 
is  as  follows:  *  5  JJ 


t,  -  pi 

ft;  da  -  vc 

yr^pz 

ftl  d$  -  V7 

0.  *  j~ 

fi  -P3 

f-Pide  -  vr 

-  V4? 
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Table  2.  Nondimensional  Buckling  Functions 

F-  0.5 

J.  5  Pi«  8.36559  P8-  4.97939  P3-  9.95878  P4-  C. 740835 
V6*  0.805893  V7-  0.43380?  VS-  0.866414  V9-  6.44004E-8 
F«  0.5 

dm  10  PI-  1.85068  P8«  3.65588  P3-  7.31056  F4-  0.41838 
V6-  0.366897  V7-  0.751816  V8-  1.50843  V9-  0.100878 
F«  0.5 

Jm  15  Pl-  1.57888  P2-  3.09786  P3-  6.19578  P4«  0.2744 
V6-  0.504807  V7-  1.08073  V8«  8.04146  V9-  0.184145 
F-  0*5 

0-  80  PI-  1.40848  P 8-  8.78485  P3-  5.56969  P4-  0.195696 
V6-  0.686739  V7-  1.86301  V8-  8.58608  V9-  0.141171 
F-  0.3 

J«  85  PI-  1.89386  P2-  8.58689  P3-  5.17857  P4-  0.143706 
V6-  0.739305  V7-  1.48808  V8-  8.97604  V9-  0.153673 
F-  0.5 

J-  30  PI-  1.81387  P8-  8*45873  P3-  4*90546  P4-  0.10608 
V6-  0.844912  V7-  1.70141  V8-  3*40881  V9-  0.168897 
F-  0.5 

Jg  35  PI-  1.15784  PS-  8.36108  P3-  4.78805  P4-  7.67278E-2 
V6-  0.945644  V 7-  1.90688  V8-  3*81363  V9-  0.169578 
F>  0.5 

jm  40  PI-  1.11987  P8-  8*89898  P3-  4.59795  P4«  5.85781E-8 
V6-  1.04307  V7-  8.10683  V8-  4.21365  V9-  0.174147 
F«  0.5 

J-  45  PI*  1.0966  P8-  8*8598  P3-  4.5196  P4-  3.1S593E*2 
V6-  1.13848  V7-  8*30343  V8-  4*60686  V9-  0*176898 
F»  0.5 

Jm  50  PI-  1.08618  P8-  8.83975  P3-  4*4795  P4-  0.01888 
V6-  1.83897  V7-  8.49889  V8-  4.9965?  V9-  0.177961 

F-  1 

J-  5  PI-  1.65491  P2-  3.53677  P3-  7.07355  F4»  1.81864 
V6»  0.143977  V7-  0.307699  V8-  0.615399  V9-  0.157699 
F*  l 

J-  10  PI-  1.42574  P2*  2.60318  P3-  5*80684  P4-  0.978366 
V6*  0.268016  V7-  0*534171  V8-  1*06834  V9-  0.248895 
F-  1 

J-  15  PI-  1*86139  P8-  2.208  P3-  4.416  P4-  0.638516 
V6-  0.377757  V7-  0.726867  V8-  1.45853  V9-  0.897846 
F-  1 

J-  80  PI-  1.13778  P2-  1*98303  P3*  3.96605  P4-  0.456095 
V6-  0.476744  V7-  0.89879  V8-  1.79758  V9-  0.3377 
F-  1 

J*  £5  PI-  1.04186  P2-  1.83657  P3-  3*67314  P4-  0.345574 
V6-  0.567386  V7-  1.05857  V8-  2.11714  V9-  0.367765 
F-  1 

Jm  30  PI-  0.965839  P2-  1.73365  P3«  3.46731  P4-  0.269123 
V6-  0.651414  V7-  1.2094  V8-  8.4186  V9-  0.391 i79 
F-  1 

J*  35  Pi*  0.904749  P2-  1.65774  P3-  3.31548  PA-  0.214054 
V6*  0.730127  V 7*  1.35362  V6-  2.70725  V9*  0.409802 
F*  1 

J*  40  PI-  0.855266  P8-  1*5999  P3-  3.19961  P4-  0.172623 
V6-  0.304535  V7-  1.49281  VS-  2.98563  V9-  0.48482 
F-  l 

J*  45  PI-  0,815082  P8-  1.55487  P3»  3.10974  P4-  0.14035 
V 6-  0.875448  V7*  1  .62809  V8-  3.25618  V9-  0.43703 
F*  1 

J-  50  Pi-  0.782553  P2-  1*51932  P3-  3*03864  P4-  0*114463 
V6*  0.94353  V7»  1.76027  V8-  3*58054  V9*  0.446988 
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r-  i.f 

dm  5  Pi,  1.60181  P8«*  8*92168  P3«  5*64383  P4«  8*93596 
06*  0,139357  V7«*  0. 854181  V8«  0*508381  V9«  0*855488 
F-  1*5 

J*  10  PI*  1.41703  P8«  8.17159  P3«  4*34317  P4*  1.50389 
V6*  0*868638  V7*  0*443108  V8*  0*886817  V9-  0*386814 
F  *  1*5 

J«  15  Pl«  1*87103  P8*  1*8573  P3*  3*7146  P4«  0*955934 
V6«  0.373818  07®  0*604693  V8®  1.80939  V9*  0*469381 
F»  1*5 

J*  80  Pl«  1*15353  P2-  1*67974  P3«  3*35948  P4»  0*669747 
V6*  0*473575  V7«  0.750831  VB®  1*50166  09*  0*587649 
Fm  1.5 

J*  85  PI*  1*05761  P2*  1.5648  P3*  3*1896  P4«  0*496506 
V6*  0.565587  V7*  0*886968  V8»  1*77394  09*  0.570845 
r*  i*5 

J*  30  PI*  0.978448  P8»  1.4843  P3»  8*9686  P4*  0*368055 
V6*  0.650718  07*  1*0161  V8*  8*0388  09*  0.604083 
F*  1.5 

.J*  35  PI*  0*918556  P2*  1.48494  P3*  8*84989  P4*  0*301336 
06*  0.730104  V7*  1*14007  VS*  8*88014  V9*  0*630343 
F»  1.5 

J*  40  Pi*  0*857384  P8-  1,37959  P3«  8,75918  P4*  0.243U6 
V6*  0.804697  V7«  1*8601  V8*  8*58019  V9*  0*651494 
F*  1.5 

J-  45  PI*  0.811017  P8*  1.34408  P3*  2*68803  P4»  0.198662 
V6*  0.875855  V7*  1*37703  V8*  2*75405  V9«  0*668778 
F*  1.5 

J*  50  PI*  0*771999  P2»  1*31556  P3-  8.63115  P4«  0.164078 
V6*  0*948419  V7«  1*49148  V8«  2*96896  V9*  0*683052 

F*  3 

*J*  5  PI*  8*13814  P8*  8*15077  P3*  4*30155  P4*  6.16254 
V6»  0*165496  V7*  0.187117  V8*  0*374235  V9*  0*536141 
F*  3 

j»  10  pi*  1*86115  P2*  1*6507  P3*  3*30141  P4«  2.80658 
V6«  0*347416  V7*  0*330729  V8-  0*661457  V9*  0*760314 
F*  3 

dm  15  PI*  1*64765  P2*  1*44892  P3*  2*89784  P4*  1*64011 
V6*  0*490761  V7*  0*456765  V8*  0*913569  V9*  0*923004 
F«  3 

dm  20  PI*  1.47617  P2*  1*33861  P3*  2*67722  P4*  1.07518 
V6*  0*619188  V7*  0*573244  V8*  1*14649  V9*  1*01654 
F*  3 

,J*  25  PI*  1*3363  P2*  1*26924  PS*  2*53848  P4*  0*754956 
V6*  0.73544?  V 7*  0*68366?  V8*  1*36733  V9-  1*08223 
F*  3 

dm  30  PI*  1*28076  P8*  1*8218?  P3-  2*44374  P4*  0. ‘>55351 
V6*  0.841653  V7*  0*78997  V8*  1*57994  V9«  1.13054 
F«  3 

dm  35  PI*  1*18431  P2»  1.1877  P3*  2*3754  P4*  0*422606 
V6*  0.939468  V7*  0*8933  V8*  1*7866  V9  *  1*16731 
F*  3 

J*  40  PIj  1*04318  P8*  1.16207  P3*  £.32414  P4*  0*330085 
06*  1.03082  V7“  0.9944  08*  1*9888  V9*  1.19602 
F*  3 

dm  45  PI*  0.974336  P2*  1*14827  P3*  2.28454  P4*  0*26303? 
V6*  1*11499  V7*  1.09378  V8*  2*18756  V9-  1.2109 

F»  3 

Jm  §6  Ptm  0-91 569 7  P 8*  1*12662  P3*  2.25324  P4*  0*818114 
V6*  1* 19465  07*  1.19179  06*  2*38359  V9*  1*23744 
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r*  t 

Jm  5  PU  3*08883  P8«  1*75614  P3*  3*51889  P4*  9*8059 
V6»  9*883456  07*  0.158785  V8*  0.305569  09*  0*853113 
F-  5 

Jm  10  PI*  8*53145  P8»  1*39831  P3«  8*79663  P4*  3*97884 
V6«  0*483698  V7*  0*874438  08*  0*548875  09*  1*1987 
r«  5 

J*  15  PI*  8*17001  P2*  1*86071  P3»  8*58143  P4*  8.13699 
V6*  0.678483  V7*  0*38418  08*  0*76884  V9-  1*38468 

f*  5 

J*  80  PI-  1*89637  P8*  1*18847  P3*  2.37695  P4«  1*314 
06*  0*837467  V7*  0.487517  V8*  0.975034  V9-  1*49893 
F*  5 

J*  85  PI*  1.68^98  P8«  1.1446  P3*  2*8898  P4*  0*876079 
V6*  0.983881  V7*  0*587098  V8-  1.17419  V9-  1*57515 

F*  5 

J*  30  PI*  1.51853  P8*  1*11554  P3*  8.23108  P4*  0.617331 
V6*  1.11547  07*  0*684149  08*  1*3683  09*  1*68886 
F*  5 

J*  35  PI*  1*37401  P8*  1*09513  P3*  8*19085  P4*  0.453034 
V6*  1*83501  07*  0*779485  06*  1.5S88S  09*  1*66887 
F-  5 

J*  40  PI*  1.85973  P8*  1.03016  P3«  8*16033  P4*  0> 34301 
06-  1*34461  V7-  0*8734  VS-  1*7466  09*  1.698i2 
F*  5 

J*  45  PI*  1*16435  P 8*  1.06884  P3«  8*13768  P4«  0*866828 
V6*  1.4459  07-  0.966389  08*  1.93878  V9«  1.78188 
F*  5 

0*  50  PI*  1.08397  P8*  1*06005  P3«  8*1801  P4-  0*210838 
06*  1*54021  V7*  1*05861  V8«  8.11783  V9*  1.73962 

F*  10 

Jm  5  PI*  4.96209  P8«  1*3953  P3«  8*7906  P4-  15*9737 
06*  0*431708  V7*  0*121391  V8*  0*842788  V9*  1*38971 
F*  10 

J*  10  PI*  3*78865  P8®  1*18467  P3«  8*36934  P4»  5*30211 
V6«  0.760793  V7«  0.884457  08-  0*448915  V9*  1*851 
F-  10 

Jm  1$  PI.  3*05189  PB*  1*11187  P3*  2.88853  P4*  8*49934 
V6«  1*08631  V 7*  0*321138  08*  0*648875  09*  2*06844 
F*  10 

Jm  SO  PI*  2.55595  P2*  1*07555  P3*  2.1511  P4-  1*39464 
V6*  1*84867  V7*  0*41471  VS-  0*82948  09*  8*18977 
F*  10 

Jm  25  PI*  2.19829  P2*  1*05514  P3*  2*11029  P4*  0*862492 
V6*  1*43993  07*  0*506508  V8-  1*01302  09*  2*26481 
F*  10 

Jm  30  PI-  1.92894  P2*  1*04289  P3*  2*08458  P4«  0*572233 
06*  1.60774  V7*  0*597187  V8-  1.19437  V9-  2*31459 
F*  10 

J-  35  PI*  1*71946  P2»  1.03364  P3*  8*06727  P4*  0*399698 
V6*  1.75734  V7*  0*687113  08*  1*37423  09*  2*34937 
F*  10 

Jm  40  PI*  1*55243  P2*  1.02752  P3-  2*05503  P4*  0.890405 
06*  1.8924  07*  0.776507  08*  1*55301  09*  2.37463 
F*  10 

J*  45  PI*  1.41664  PC*  1.02302  P3»  2.04605  P4*  0.81766: 
06*  2.01564  07-  0*86551  08*  1.^3102  09*  2.39357 
F*  10 

J*  50  PI*  3.3045  P2*  1.01963  P3*  8.03926  P4«  0*167342 
06*  8*12914  07*  0.95481S  02*  1*90844  09*  2*40813 
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